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Abstract 

In the rough phase, the width of interfaces separating different 
phases of statistical systems increases logarithmically with the sys- 
tem size. This phenomenon is commonly described in terms of the 
capillary wave model, which deals with fluctuating, inflnitely thin 
membranes, requiring ad hoc cut-offs in momentum space. We inves- 
tigate the interface roughening from first principles in the framework 
of the Landau- Ginzburg model, that is renormalized field theory, in 
the one-loop approximation. The interface profile and width are 
calculated analytically, resulting in finite expressions with definite 
coefficients. They are valid in the scaling region and depend on the 
known renormalized coupling constant. 
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1 Introduction 

Interface roughening is a phenomenon which has attracted interest of ex- 
perimental and theoretical investigators, see e.g. [Ij EJ [3l HI [3 O [7] , since 
its discovery [8 . It is displayed by interfaces, separating different coexist- 
ing phases or substances of a system of statistical physics, in a range of 
temperatures Tr < T < Tc between the roughening temperature Tr and 
the critical temperature Tc- Roughening manifests itself in a characteristic 
dependence of the interface width on the system size. For an interface of 
diameter L the width increases logarithmically with L in the rough phase, 
whereas it remains constant of the order of the correlation length ^ for 
temperatures below Tji. 

This effect is commonly described theoretically in terms of the capillary 
wave model or drumhead model [8J. In this model the interface is repre- 
sented in an idealized way by an infinitely thin fluctuating membrane, so 
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that the instantaneous microscopic interface profile is a sharp step function 
between the two phases. Nevertheless, in the thermal average the capillary 
wave fluctuations produce a continuous density profile with a finite width 
w, which can be shown to be given by an integral over all wave-numbers of 
the fluctuations, which is essentially of the form 
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where D is the number of dimensions of space and a is the interface tension 
A natural lower limit on the wave numbers is given by the system size, 

const. 



L 



(2) 



In order to avoid the divergence of the integral, an upper cut-off fcmax has 
to be introduced. As there should be no waves with wavelength smaller 
than the intrinsic width of the physical interface, the upper cut-off is taken 
to be of the order of the inverse correlation length. In the case of Z? = 3, 
considered here, one obtains 

w'^-^\n^ (3) 
ivrcr 

with an unknown constant c. The logarithmic increase with L is due to 
the contribution of capillary waves with long wavelengths near the system 
size L. 

Complementary to the capillary wave model is the mean field descrip- 
tion of interfaces. In mean field theory and its field theoretic refinements, 
interfaces possess an intrinsic continuous profile with a well-defined width, 
which is proportional to the bulk correlation length and does not depend 
on the system size. 

Mean field and capillary wave theory can be combined in the "convo- 
lution approximation" [9l [TOl . In this picture the intrinsic profile describes 
the interface on a microscopic scale of the order of the correlation length, 
while capillary wave theory describes the macroscopic interface fluctuations 
of wavelengths much larger than the correlation length. The intrinsic pro- 
file is thus centered around a two-dimensional surface subject to capillary 
wave fluctuations. In the convolution approximation the square of the re- 
sulting total interface width is obtained as a sum of the intrinsic part and 
the capillary wave contribution, 

^2 = Cie^ + -^lnA. (4) 

The description of rough interfaces by means of the capillary wave model 
and the convolution approximation is unsatisfactory for different reasons. 
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First of all, it has so far not been possible to define the concept of an 
intrinsic interface profile and width unambiguously outside a given theory. 
In experiments or Monte Carlo simulations of systems with interfaces, the 
observed interface profile and width are the total ones, including the effects 
of the intrinsic structure as well as of the capillary waves, and there is no 
clear way to separate the intrinsic structure from the effects of capillary 
waves. Secondly, the models sketched above contain ad hoc constants, 
whose numerical values are arbitrary and cannot be fixed unambiguously 
within the models. 

In this article we investigate the profile and width of rough interfaces 
in a coherent approach from first principles. Statistical systems with co- 
existing phases, separated by interfaces, are described in the framework of 
the field theoretic version of the Landau-Ginzburg model, including fluctu- 
ations on all length scales. No cut-off on wave-numbers is introduced. For 
explicit calculations we employ the one-loop approximation. It should be 
noted, however, that an extension to arbitrary higher loop orders is possible 
in principle. The interface profile, resulting from the calculation, shows the 
expected logarithmic broadening with the system size L. We obtain ana- 
lytical results for the numerical coefficients, which are fixed unambiguously 
in this approach. 

Interfaces have been studied before in the framework of field theory by 
other authors. In [TTl [T^] the profile is calculated to first order in the e- 
expansion, where D = 4 — e and an extrapolation to e = 1 is necessary. The 
e-expansion is an expansion around the four-dimensional case. As can be 
seen from Eq. ([T]), in four dimensions the contribution of long- wavelength 
modes converges and no roughening is present. This has the consequence 
that within the e-expansion, even after extrapolation to D = 3 dimensions, 
roughening effects do not show up, as is well known. The calculation of 
[TT] is extended to include the effects of an external field in [T5' . 

Our calculations are performed in 13 = 3 physical dimensions in con- 
trast to the e-expansion. The three-dimensional approach is based on a 
systematic expansion in a dimensionless coupling [14|, I15j . Ultraviolet di- 
vergences are treated by dimensional regularization (D = 3 — e), which 
does not vitiate the fact that the results for physical quantities strictly re- 
fer to D = 3 dimensions. This is also seen explicitly by the fact that the 
calculation reveals the typical roughening effects. Renormalization of the 
three-dimensional field theory is performed in the scheme used in [16j to 
two- loop order, employing the results of [TTl fTS] . 

A three-dimensional study has previously been done in [1^ . where the 
interface profile is considered in Z) = 3 dimensions at one-loop order in the 
presence of an external gravitational field. A functional form of the profile 
is given, including capillary wave effects. The dependence on the system 
size is, however, not considered. We shall compare our results with the 
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ones of [H] below. 



2 Interfaces in field theory 

In the framework of field theory, the system under consideration, possessing 
interfaces, is described by an order parameter field representing the 
difference between the concentrations of the two coexisting phases. The 
physics of the system is governed by the Landau-Ginzburg Hamiltonian [20] 

= Jd^xn{(l){x)) (5) 

with the Hamiltonian density 

Hi(l))^^d^^d^^ + Vo{(t>). (6) 

In the situation with interfaces the potential is of the double-well type, 

VoW = ~ vl)\ (7) 

Mean field theory amounts to the classical approximation where fluctu- 
ations are neglected. The minima of the potential then correspond to the 
two homogeneous phases. The mean field correlation length is defined 
through the second moment of the correlation function in the mean field 
approximation. It is given by the second derivative of the potential in its 
minima: 

Co = (K'(«o))"' = A- (8) 



With the bare mass mo, defined by 



mo = 7", (9) 
the Hamiltonian density can be written as 

= - !^<^2 + ^ M . (10) 

The simplest description of interfaces is also based on mean field theory 
|21j . In this approximation the interface profile is given by minimization of 
the Hamiltonian H with boundary conditions appropriate for an interface. 
The corresponding field equation 
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leads to the differential equation 



A(/.-l^o'W=0. (12) 

If we choose the interface to be perpendicular to the z-axis, we find the 
typical hyperbolic tangent profile [22] 

4^°)(z)=«otanh(^i^). (13) 

Its width is proportional to the mean field correlation length ^q. The 
parameter zq specifies the location of the interface. 

Essential for a field theoretic treatment, as being considered in this 
article, are corrections to mean field theory coming from fluctuations of 
the order parameter field. They can be calculated systematically in renor- 
malized perturbation theory. The fluctuations result in different modi- 
fications of the mean field result, as will be considered in detail below. 
First of all, higher order corrections change the form of the profile from 
the tanh-function to a different function. Secondly, renornialization of the 
parameters vq and becomes necessary and, as a result, the mean field 
correlation length is replaced by the physical correlation length ^, which 
diverges near the critical point with a characteristic exponent v. Finally, 
long-wavelength fluctuations lead to the roughening phenomenon, which 
implies a broadening of the interface, such that its width depends logarith- 
mically on the system size and diverges in the limit of an infinite system. 

The partition function for the system with an interface can be written 
as a functional integral of the form 

Z = JvipcM-H[^o + V']), (14) 

where (j^o^z) is a classical interface solution as given above and ip{x) denotes 
the fiuctuations around it. The Hamiltonian density, expressed in terms of 
(fi, reads 

■Hicj^o + ^)= n{M + lip{x)K^ix) + + ^^\^) • (15) 

Here the operator K is given by 

if--A-^ + f^^(-), (16) 

where A is the Laplacean. 

In the loop expansion the quadratic terms in Ti. are treated by means 
of Gaussian functional integrals, and the higher order terms are taken into 
account by Taylor expansions. 
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The spectrum of K is known analytically [23]. We have to employ it 
for our calculation and give details below. At this point we would like to 
draw the attention to the fact that K has a single zero mode 

Kil){x)^Q. (17) 

The zero mode of the fluctuation operator is directly related to translations 
of the interface, as parameterized by the parameter zq. For every value of 
this parameter, the hmction ^o^""* is a solution of the classical field equation. 
This implies that 

M.)^^ (18) 

is a zero mode of K. 

The existence of a zero mode requires to treat the corresponding fluctua- 
tions, which are proportional to ■0, separately from the remaining Gaussian 
integrals in the functional integrals. This is done by the method of collec- 
tive coordinates [21] • The collective coordinate in question is zq. In the 
Gaussian integral it is set to an arbitrary value, which we choose to be 
zq = 0, and the fluctuations are restricted to the space A/j_ of functions 
orthogonal to the zero mode V'- 

dPxip{x)i}}{x) ={). (19) 

When expectation values in the presence of an interface are calculated, 
integration over zq would imply averaging over all translations of the in- 
terface, leading to translationally invariant results. In case of the interface 
profile, however, this is obviously not appropriate, since one is interested 
in the profile function relative to the position of the interface. Therefore 
integration over zq ha-s to be omitted, leaving us with Gaussian integrals 
over AG. . So the interface profile is given by 

</>e(a;) = <^o(x)-f <^/(x) (20) 

with 



<j,f{x) = {^{x)) = Vip^ix) exp(-iJ[0o + V]) ■ (21) 



3 The profile equation 

For functional integrals over the fluctuation field ip G J\f±_ Fcynman rules 
can be set up analogously to the usual case. The propagator and vertices 
can be read off the Hamiltonian . The propagator is the inverse of the 
fluctuation operator restricted to A/j_ : 

K'^K\^^. (22) 
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There are three-point and four-point vertices, given by 




= ~9o- 



The fluctuation part the interface profile gets contributions from all 
orders of the loop expansion: 



= (l)i{x) + 4>2{x) 



(23) 



In the one-loop approximation, which we employ, the Feynman diagram 
contributing to the profile function 




leads to 



(j)i{x) 



90 
' 2 



S'x' K'-^(x., x') K'-^{x', x')Mx') 



(24) 



Here the kernel of the inverse operator K'^^ enters. It would be possible 
to calculate 0i from this expression. It is, however, more convenient to 
obtain it as a solution of a differential equation. Acting with the operator 
K on Eq. (|24|) . we obtain the profile equation 



90 



^{x,x)(j)o{x) = 0. 



(25) 



In order to solve this equation we need the explicit form of K'^^{x, x), 
which is discussed below. 

An alternative derivation of the profile equation is based on the so-called 
effective action r[$], which is a functional of a field ^{x). r[$] is obtained 
by Legendre transformation from the free energy in the presence of a non- 
constant external field. For constant $ the effective action reduces to the 
Gibbs potential. For a definition and discussion see e.g. [25]. Calculating 
r[<I'] in the one- loop approximation and finding the interface profile 4>{x) 
as a stationary point of F, 

0, (26) 



6<^{x) 



again leads to Eq. 
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4 Solution of the profile equation 

The inverse of the fluctuation operator K' at coinciding arguments, which 
enters the profile equation, can be obtained by means of the spectral rep- 
resentation. K is the sum of the negative two-dimensional Laplacean and 
a one-dimensional Schrodinger operator 

K = ^A'^^'>+K, (27) 

where 

K = -d! + ml-^-^sech'('^z). (28) 



2 V 2 

The negative Laplacean on the L x L square has eigenvalues 

with k = — n, n G Z^, (29) 
and corresponding eigenfunctions 

^^i(f) = L-i e'*"^ xe[0,L]\ (30) 

The spectrum of K is known exactly [23 . It consists of two discrete eigen- 
values 



Lo^'> = 0, Mz) = \/ ^ sech^ [-Y^ ) , (31) 



a;W = Jm^, (,) = tanh (l^zj sech (1^. ) , (32) 



/3mo ^ , /Too \ , /"rno 

^.-..(j, v^iy-y y TS\n h ( 7, I fipr h ( 

and a continuum 



ojp — rriQ + p with p G M, (33) 



2p H TOq tanh I — z I + diTOoptanh I ;; 

2 2 V 2 / V 2 

(34) 

with the normalization factor 

AAp = (27r(4p4 + 5m^p2+m4))-i (35) 

The spectrum of K is thus given by 
47r2 

AfScj = -p-n +u!, 'i'iiujix) = (pn{x)ipuj{z), (36) 

where w runs through the eigenvalues of K. The zero mode, discussed 
above, is represented by '^^^q- 
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In terms of the spectrum we write 

K'-\x,x)=^Mx)rx{x)\. 

A 

Inserting the exphcit expressions, we obtain 



(37) 



K'-\x, x) = Co + (Ci -C2 + Ci) sech" (^^) + {C2 + C3) sech^ (J^z) , 

(38) 

where the coefficients Ci are 



Co 
C2 



n 

y— 



27r 

3mo 



47r2n2 + (mg+p2)L2 



3mo 



E 



47r2n2 ' 



4 4^ 47r2n2 + |m^L2 ' 



C3 = -3m^ 



niQ +p2 



n 

C, = JdpM^ ^ W + (m2+p2)x2 



1 



(39) 

(40) 

(41) 
(42) 
(43) 



These expressions are divergent and have to be regularized, as discussed 
below. 

With the explicit form of K'~^{x, x) at hand, the solution of the profile 
equation is found as 



9oVo 
2ml 



Cotanh (^2:) 
(Ci-C2 + C4)tanh(^^) 



(Co + C2 + C3)^Z 



scch2 ( — -2 
\ 2 



(44) 



Written in this way, the expression for the profile contains the divergent 
coefficients Cj as well as the bare parameters go, and vq. In order to ar- 
rive at a finite expression in terms of physical parameters, renormalization 
has to be performed. 

The divergences have to be treated in some regularization scheme. We 
choose to employ dimensional regularization in Z) = 3 — e dimensions. 
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It should be noted that this does not amount to an e-expansion, since 
after renormahzation e is sent to zero, whereas in the e-expansion one 
has D = 4 — e and the results have to be extrapolated to e = 1. So 
our use of dimensional regularization does not vitiate the fact that the 
results for physical quantities strictly refer to D = 3 dimensions. Using 
other regularization schemes, like Pauli-Villars, would lead to the same 
final results. 

We adopt the renormahzation scheme used in jTFj to one-loop order. 
The renormalized mass mn = 1/^ is equal to the inverse correlation length 
5, which in turn is defined through the second moment of the correlation 
function. The field and its expectation value v are renormalized according 
to ^ ^ 

(I>r{x) = ^= (l>{x) , vii = —=v, (45) 



'Zr, V Zr 

where Zr is the usual field renormahzation factor. The renormalized cou- 
pling is specified as in [23] through 

9R--f- (46) 

"r 

In addition we define a dimensionless renormalized coupling according to 

UR = • (47) 

^R 

Employing the relations given in [171 118] , the bare quantities mg and go 
are expressed in terms of their renormalized counterparts. 

The coefficients C'i are evaluated in the same scheme. Leaving out the 
lengthy details, we quote the results 

Co = (48) 



C2 + C, = ^\n3 (49) 
Ci-C2 + Ci = ^(-a + In(moi)) (50) 



-Inl^;^ 1-7- 1.832, (51) 



with 

n/ — In I 

where 7 « 0.577 is Euler's constant. 

The coefficients Ci contain additional terms decaying exponentially fast 
with L, which we neglect here. 

Inserting everything into the expression for the interface profile yields 
the renormalized interface profile 4>r{z) depending on the parameters mR 
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and ur. In this expression the divergences are cancelled, as they should. 
Expanding consistently in powers of ur up to the first order, we obtain 

(I)r{z) = Wflj tanh(^z) 

H (a — ln(m/;L)j tanh( z) sech ( z) 

167r 2 2 

which is the central result of this article. Asymptotically, for \z\ oo, the 
profile approaches the bulk expectation value ±vr, corresponding to the 
pure phases of the system, as it should be. Note that the profile depends 
logarithmically on the system size L, revealing the effect of capillary wave 
fiuctuations. It is this term, depending on m/jL, which represents the 
deviation from the Fisk-Widom [27j scaling form (f>(mRz). 

In order to illustrate the characteristics of the interface profile, we have 
to specify a numerical value for the dimensionless renormalized coupling 
ur. In the vicinity of the critical point the coupling varies only slowly 
and is close to the universal fixed point value = 14.3(1), see [55] for 
a discussion of numerical and ficld-thcorctical estimates. Therefore we 
take Ur — 14.3 in the plot. The interface profile according to Eq. (I52|) is 
displayed in Fig. [T] for different values of mRL. 




Figure 1: The renormalized interface profile for difi'erent system sizes at a cou- 
pling of Ur = 14.3. 
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For L larger than 200 ^ the profile is no longer a monotonic function of 
z. In this region the one-loop contribution approaches values of the size of 
the leading order term and the one-loop approximation reaches its limit of 
validity. 

A three-dimensional calculation of the interface profile by field theoret- 
ical methods has been performed previously by Jasnow and Rudnick [19] 
in the one-loop approximation. They do not consider systems of a finite 
extent L, but use an external gravitational field in order to control capillary 
wave fluctuations. 

Their result for the interface profile has the functional form of Eq. (|52p , 
with coefficients that are given numerically. In place of our term 



^ln(m7jL)tanh(^z)sech^(^z) (53) 



they get 



C3ln(/i)tanh(^z)sech2(^z), (54) 

where h is the external field, and C3 — 0.109975. As this term originates 
from long- wavelength fiuctuations, its coefficient can be expected to cor- 
respond to ours. For a comparison one has to take into account that h 
corresponds to 1 . Also, they use an estimate of the fixed point value 
of the coupling which is slightly different from ours. In view of this, the 
numerical coefficient C3 is in rough agreement with our result. 

The coefficients of the other terms are different from ours. We do not 
know, whether there is reason to expect them to be the same. 



5 Interface width 

There are various ways to define the width w of an interface, see e.g. refs. [9] , 
[5] and [55]. A suitable choice is 

= (z^) = Jdzz^p{z), (55) 

where the weight p{z) is taken to be proportional to the square of the 
gradient of the profile, 

p{z) cx {dMz))^ (56) 



dzp{z)^l. (57) 



and to be normalized: 



In the evaluation of w'^ it should be observed that the occurring terms have 
to be expanded consistently in powers of the coupling. For example, in the 
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one-loop approximation, the square of the gradient dz4> — dz4>o + dz4>i is 
to be taken as 



(58) 



With the interface profile given above, we obtain 



w 



m|j 207rTOjj, 



with 



12 



(TT^-e) ln3 



13 
12 



(59) 



(60) 



The interface width grows logarithmically with the system size. So 
the field theoretic calculation in the one-loop approximation confirms the 
prediction of capillary wave theory. Our result, however, does not rely 
on the capillary wave approximation, but comes from taking into account 
fluctuations of the density profile on all scales. Moreover, the numerical 
coefficients are fixed unambiguously and do not depend on ad hoc cut-offs. 

For a direct comparison with the convolution approximation another 
definition of the interface width is more convenient, namely choosing 



p{z) = {2vRy' dz(t>R{z) , 



(61) 



which is meaningful as long as the profile function is monotonic. For this 
choice, in the convolution approximation the squared width w'^ of the in- 
terface equals the sum of the intrinsic and the capillary wave contributions, 
see Eq. 

With this definition of the interface width, our result reads 



Ur 



w — — n- -I- 5-ln(mfli) 



with 



ur r. 2A ^ 13 



1.249. 



(62) 
(63) 



By noting that for the interface tension we have [TB] 

1 Ufi 



a 2to|j 



-OK), 



(64) 



we see that the L-dependent term is in agreement with the prediction 
from the capillary wave model. 

The L-independent term a contains the classical mean field value 



ao 



(65) 
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see |32| . plus corrections, which are undetermined in the convohition ap- 
proximation. 

For both choices of the weight function ^(z), the one- loop approximation 
ceases to be valid, if L gets so large that p{z) becomes negative. This 
happens for L « 200^, which coincides with the value, where the profile 
begins to be non-monotonic. 

Interfaces have been investigated in the three-dimensional Ising model 
by means of Monte Carlo calculations in [SUl [13 [SB [13 • To observe rough- 
ening in Monte Carlo is delicate, nevertheless one can obtain estimates for 
the offset a from their data, which amount to a = 2.68 [30], a = 0.76 [29], 
a = 3.44 [31] and a = 0.08 [32]. In view of the spread of these numbers, 
and in view of the fact that higher-loop contributions will change our es- 
timate, we can only notice that the order of magnitude is compatible with 
our result. 

6 Conclusion 

Field theory, in the form of the Landau- Ginzburg model, including thermal 
fluctuations on all length scales, allows to determine the interface profile 
and interface width for models in the Ising universality class in the critical 
region. For a system possessing a square interface of size L x L we derived 
the conditional equation for the interface profile in the one-loop approxi- 
mation. We obtained its solution in analytical form. When it is expressed 
in terms of physical, renormalized parameters, no divergences occur, and 
there is no need to introduce ad hoc cut-offs as in the capillary wave model. 

The solution displays the characteristics of roughening by depending 
logarithmically on the size L. The interface width grows logarithmically 
with increasing system size. The coeflScient of the logarithmic term is in 
agreement with the universal part of the capillary wave model, and the 
constant term is consistent with results from Monte Carlo simulations of 
the Ising model. 
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